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IXTHRODUC TI ON 


Today the subcritical assembly is a convenient tool of 
the nucletr seientiet and engineer. It hes several inherent 
edvantages which weke it particulerly useful. It is much 
Cheaper to sssemble end operete then is a full seale reactor. 
The amounts of wmeterials used ere much less. It is inherently 
sele, no extensive shielding need be provided, snd work on 
be carried on with none of the physics1 restrictions neces- 
Bary Neer a resctor. 

A netursal uranium, graphite moderated suberiticel sseen- 
biy hes been constructed recently at Iowa Stete College. Thie 
assembly will be useful in sat lesst two caprecities. First, 

4% will provide extensive instructional opportunities, and 
the chence to compare clossrocom calculistione with eetual 
experimental results. Secondly, it will provide a reeeerch 
faciiity for numerous thesis investigstions by grecuete 
stucents. 

Before other studies eculd be underteken, it was neces- 
Bary to develop eefe and efficient procedures e¢ well eas to 
Getersine the besic operating chsracteriatics of this assen- 
bly. It wae with this .urpose in mind that the present thesis 


wes undertezen. 
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REVIEW OF Tat LITERATURE 


The histories] role pleyed by tre subcritics] aegsembly 
in achieving the firet self sustaining nucleser chain resection 
hae been given by Smyth (1). It is summerized below. 

On December 2, 1942, mankind initisted his first self 
sustaining nuclerr crrin resection at the tinivereity of Chi- 
¢ago (1, p- 7). This event merked the culminstion of months 
of concentrated effort by meny scientiets and engineers. It 
@eculd not heave occurred as it did without the suberiticsl1 
assembly leeding the way. 

By June 1940 it hed been Getermined thet the best pos- 
Sibility for & self sausteining chein revetion lay in the 
thermal neutron Pission of yess | using a heterogeneous wix- 
ture of graphite and naturel uranium (1, p-. 38). It was 


375 5 
known that in the fiesion of vy”? 


more than 2 fast neutrons 
were produced, on the aversge, for every neutron consumed 

lun Cission. however, moat of theece aepeare neutrons were 
absorbed by other meterisis in slowing down, or escaped alto- 
gether, and thugs were lost before they in turn were able te 
react with other uranium nuclei to produce fission. The pvrob- 
lem wes to errange the flesioneble moterial end moderetor in 
such a manner that enough of the neutrons were conserved to 
ayuasteain a shein resction. 


Two mein avenues of approach to this problem were por- 


Sitle.- An elaborate series of investigetions into the nuclear 
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properties of verious meterirels wea initisted, end theoretical 
celeulations were begun besed on these deta. towever, since 
speed was the controlling fretor et thet time, the major 
effort wes directed into the empirical approach, and here 

the suberiticr)] nasenbly cecauwe the primery mesns of investi- 
geting the croblen. 

The very emell emounte of meterinis available, enc the 
extreme costs involved, made full senle experimente ispos- 
sible. The method used then was to eonatruct a plle of ebout 
oue fourth to one third the size believed necessary for 
eriticality. 

From the results of experiments conducted with such a 
Subcriticel essembly it wes possible to predict the chrrecter- 
istics of a similar full seale reeaetor. 

sinee such ~@ euberiticel essembly could not sustsin a 
ehain resection an exterrel source of neutrons wre pleceé near 
the bottom of the atructure. Tre flux distribution in sueh 
en eseembly then decreased exponentially with distenee from 
the external source. This feet led to the subcriticr] assem- 
Bly being celled the “exponentirl pile®. | 

Doctors Fermi and Szilard hed euggested the use of 
eraphite as a moderator for e chain resction. Further, it 
was they who developed tre lattice structure, in which lumps 
ef urenium would Ge placed at reguler intervels in a matrix 
of grephite (1, p- 23). Thia lattice structure bed definite 


Avantages over e uniform sixture of uranium and moderator. 
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In July 1941 the first subcritical essembly was set up 
at Columbia University. It wee a graphite cube & ft. on a 
side, containing approxiustely 7 tone of urenium oxide dis- 
trituted in lumps throughout the graphite. An eternal 
radium-beryllium neutron source weg placed necr the bottom 
of the aseembly. Similer structures were set up in September 
end October of 1041. The multiplicetion factor, k, and the 
infinite multiplicetion factor, k.,, were computed in exch 
ease. For the third essembly, k., rerohed 0.87 (1, p. 48). 
Thwe problez remaining was to get k,, to remch or surpass 
1-00 Gy welne ~urer meterisls end possibly by using different 
lattice arrangements. 

#ithin the following yerr a series of experiments was 
eonducted using suicriticel agsemblies. In ereh experiment 
the uraniupe-graphite ratic was varied, oxides of imroved 
purities were used, lattice spacing was varied, and vericus 
Sizes and shapes of fuel lumps were tried. The ninth pile 
of this serles was constructed in July 194&. The graphite 
used in tris cese wae sufficiently free of Impurities thet 
it aveorbed 20 per gent lesa neutrons than the best previously 
avelileble commercial graphite. Calenlationse bsaged on this 
essexubdly showed that for s reector of infinite size, k = 
1-007. An infinite multipliceticon feector greater than 1.00 
hea been reached (1, p-. 69), and tre suberitienl sasembly 
head teen the mesns of rehieving this goel. 


Soinee 1942 the suberiticnal aseseably has been used 
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extensively in developing new reector types end in studying 
aesocisted problems. At tie Hanford Atomic Pro¢ucts %perstion 
over 100 suberitics] experiments have been conducted (?, p. 
303). These experimenta have had as their objective the find- 
ing of ae better water-cooled granhite lattice for plutonium 
production and providing improved knowledge of lettice theory. 

Argonne National Leboretory hse been conducting a series 
of fact exponentiel experinents (2, p- 342}. This series of 
esperimente wes cesi,ned to obtesin funcdementsl information 
about the physics of dilute fant reactors. In connection 
with the Cehool of Nuclesr Science and Engineering, Argonne 
Rational Laboratory also operrtes o ywraphite moderated netural 
uranium subcriticel assembly es well ss a heavy water moder-~- 
sted natural urenium asseably. 

in omer to ald in the design of liquid metal cooled 
reactors, North American Avietion has conducted a series of 
Subcriticel experiments with a simulsted sodium coolent (2, 
p- 309). these experiments included both natural and en- 
riched urenium fuel elements in « graphite moderator. Korth 
Americen Aviation hes aleo conducted a series of suberiticnl 
experiments to investigate the besiec physics of a heavy weter 
moderated lattice using vsrious fuel enrichments (2, p- 268). 

Beginning in 1944 and 1945 a series of subcritical 
experiments was ccnducted at the Oak Ridge National Lebora-~ 
tory. The purpose of these experiments was to determine 


Values of buckliny end other reloted cheracteriatics for 
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verious Giameters of naturel uranium fuel rods in combination 
with a series of Gifferent volume retios of light water 
moderator to uranium (2, p- 185). At the present tise the 
Oak Widge School of ANeector Teehnology hes a graphite mocer- 
ated neturel urenium suberitical assembly for use primerily 
as an inatructionel tool (3). 

A unicue applicrtion of the suberitice] sasembly was 
mede curing the construction of the Brookhaven rescter. 
During, a two-month period the pertially completed reeetor 
wee used for s series of subcritical me»surements. The re- 
@ulte of these measurements were used to determine critical 
buckling end other Lattice constente for the full serie 
reactor (2, p- 305). Later a series of suberiticsl! experi- 
ments wee conducted sot Brookheven to investigate the char~ 
acteristics of lignt water modersted rerectors with slightly 


enriched fuel (2, p. 184). 
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TAMORY GF THY SORCRITICAL ARDEMBLY 


%@ flux Gistripution in s suberitical sesembly does not 
Gatisty the wave @yjuation for na eritics. reretor. soavever, 
if the eseembly is quite large the thermal meutron flux dis- 
tribution at a distenee from beunderies anc from the scuree 
Gai ce represented quite well ty the equetian 
VP 8 P= (1) (4, aan 
BS is the materiel buckling for the systes under coneldern- 
tion. Although this equation is valid only for ec homorenecus 
gyeten, Lt may be used for a haterorpeneous system without 
gericus error. There will be loenilszed irregulsraties, but 
Zuuatien 1 will give the overall neutron Cistribution. 
weteriel buckling is determined from sersuremente of the 
therfel neutron flux distribution throushunt the suberi tier) 
aesenibay. Tre gritieal size of ¢ renetor of the anne compor)- 
tion end geometric configuraticn een then be deterrined by 
@queting gevmetrige] tuekling to the seteriel buckling ¢eter- 
mlned ebove. 
for the usual rectangular suberitics] assemcly "quation 


i i8 expressed as 


dP, 22°F, 92, ye 2 © (2) 
Ox” dy dn 


ere consijered to Le sepereble verlebles, 
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ZquAation 4 can oe written 
ata Be yt a nk (8) 
if the origin of the coordinete system le ssaumec te de 
et tree center of the side of the sasembly neerent the external 
Reulron soureée, tee boumdsry conmlitions gen ve exnressed ae; 
it- the thermal neutron flux due te the extermiu)] source 
of feet neutrons equala s conatent, Po, “ver the 


eatire plane 2 = U, em 


a+ Ttherme] neutron flux p O et xs + &/2, y = = b/8, 


ena g2 = @. (a, b, end © ore overall dimensicna o 
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the sssembly, Lucluding extrasoletion distence.) 
Fer amy given Value of 2, te flux is es mexieam at » 2 4, y 2 
O, end ia aymmetricelly dietributed. It thea follows thet 


Xm @ Ap, coe (21x) (2) 


Yn = Cn eos (ST) (10) 
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£ TI 2 
o® 7 (2) (11) 
and Bo 3 (ary (12) 


Substituting in Equation & leads to the form 


x ee (amyh + Bt) * (13) 


The solution for Z is of the form 2(,) = F sinh d(e - 2), 
and the general solution may be expressed ss 
Zon * Fan 81h Sanle ~- 2) (14) 
Combining Faustionsa G, 10 end 14, the general solution of 


the wave equation becomes 





PB x.y,2) a 7 hon cos( TS) cos( BEY) sinh Vun(e - 2) 
(15) 


Sd 


§ 
where Ay, Cp, ond Fen are combined into A,_- 


If the external neutron source is at leest two slowing 
Gown lengthe below the z = © plane, then it can be assumed 
that virtually all neutrons have been thermalized ty the 
tice they enter the subcritical assemtly. Under thoee con- 
Gitions the only term of Equation 15 thet ia of any signifi- 
Gence is the fundamental mode (m= 1, n = 1). Pauetion 1é 


then reduces to 


p. Ay Cos HS cos ro sinh Sle > 2, ) (18) 
P. ..) = F sinh 3(e - «) = g* le Blc-z) _ en B(o-2} (17) 


Tris can be further reduced to 
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Po) = F'' - be fo~ ¥ ™ e72 te 0 ef (18) 


For given physical dimensions, e FF 45 constant, end 
Equation 18 beesnes 
PD») _ cq! e~ 2 E a o-2 F (on) (19) 
If the verticel dimension, ¢c, is very lerge, end if 2 1s 
not esllowed to eyproseh c, then [; - e-2 ¥(en2)] approaches 
unity. Equetivn 19 simplifies to 
Pix) = 0 om F (20) 
Equation 20 indicates thet tre thermel neutron flux distribu- 
tion Gecreases exponentially wits increasing z. Thue e plot 
of diatence « versus the naturel logarithm of the flux (or 
of a quentity proportional to flux) should te a straight 
line of slope ¥. Neer the top of the essembly the points 
Will Geviste from a streight Line due to neglecting the end 
ecrrection tern fa - e~2 Fez) in Tquetion &0. 
The slope, Sin ie reslly Sai: and Equetion 15 can be 


written 
va re & Q 
(It) + () sy” RY (21) 


This eguation can oe solved for materiel buekling, all 

Qther quantities now being known. Gy setting thie value of 
Material buckling equal to geometrical buckling, the critical 
@ize of the full ecele counterpart of t>e subcritical asses- 


bly can be ceternined. 
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DESCRIPTION OF APPARATU® AND EAPERIMENTAL PROCEDURE 


Figure 1 is an oversll view of the suberiticn] assembly. 
Figure < shows importent dimensiones and the errengement end 
numbering of the foll slots. The rxasexbiy consists of 
cylincers of reactor grade graphite 60 in. long ané 7 in. in 
Giameter. These cylinders hed the sides cut flirt to form 
@gurres 6 in. ecross with rounded corners. There sre steeked 
in a equare lettice, ten rows wide end fourteen rows high. 
The top five rows sre only & in. high by 6 in. wide, making 
the overall dimensions of the assembly 60 by 60 by 79 in. 
nigh - 

This plle is get on se wooden bese, anc an externel 
meutrcn ecurce is centered under it. Vater filled aluminunc 
tanks entirely surround the external neutron source, provid~ 
ing adequste shielding egsinst both gamma end neutron redia- 
tion. 

The rounded corners between the graphite blocks presen: 
117 holes spproximstely 11/2 in. across for insertion of 
fuel elements. The ospeeing of these holes allows a mitimun 
fuel lattice of 6 in. when all holes are filled. By remov- 
ing alternete fuel elements, or entire rows of fuel elements, 
the lattice spacing may be varied from 6 in. te 8.5 in. 
(6.484 in.), 12 in. or larger if desired. 

Thirty-one small horizontel slots pierce the rssesbly 


from the erst face (Figure 2). These slots provide for the 
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Figure 2. Subcritical assembly - dimensions and 
foil holder errangement 
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irradiation of indium folls to use in determining the neutron 
flux at Verious points throughout the pile. One vertical hole 
ie available for possible eimuletion of control rod effects. 
All vertical sides and the top of the aseembly sre covered 
with 0.010 in. sheet cadmium to present essentislly a “black 
boundery* to neutrons. 

The indium foils ueed in this experiment sare 1.0" in. by 
1.50 in. by 0-005 in. thick. They ere mounted on sluminum 
planchets with DuPont "Duco" cement, eseh with the weight and 
foil number inscribed on the back. The average weight of 
the 160 foils in use is 0.5945 graeme. Aluminum foil holders 
ellow placing of the foils in the previously “cescribed slots 
at 5 in. center to center spacing scross the sacerbly. For 
meesuring fest neutron flux indium folls esn be wrapned in 
O.010 in- eadmiun sneet. 

Foil counting epperatus consists of s suliteble decimel 
scalar and either a mica end window Gelger-Nueller counter or 
a thin glass wall Gelger-Nueller counter. 

Tne use of indium folls to messure neutron fluxes makes 
use of the so-called activetion method. In this method the 
indium ia exposed in the neutron flux. Indium hes a high 
cross section for thermal neutrons and becomes radloretive, 
emitting e beta particle with a 54 min. half life. The 
neutrons are quite difficult to mersure directly, but the 


beta activity of the indlum is proportional to the neutron 
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activity, and bete particles can te detected and counted 
quite readily. 

For eech experimental run indium folls were placed in 
the essembly and irredisted for a minimum of 6.4 hr- or ealx 
half lives; in most cases the irredirstion time wes much longer 
tran this. This procedure insured an activity of et lerst 
98.4 per cent of maximum of the 54 min. activity. Foils were 
then retoved individuelly or in groupe of not more then 
@ight, depending upon the perticuler experiment. The foils 
were counted in «a nesrby room heving * much lower backpround 
than the room containing the suboriticel essembly. Fach foil 
wags counted for 5 min. Then these ecunts, less background, 
were corrected back to the tine of removel from the sasrexbly 
to give net sctivity in counts per minute- In esc: instence 
in which more than one foll was removed for counting «st ea 
given tize the foil expected to heve the leest activity was 
counted firat. This helped to minimize distortion of low 
counting retes by large corrections to time of rewtovel. 

Ho foll wae counted sooner than 3 min. after resoval fron 
toe assembly. This allowed a 15 sec. indium eetivity to 
decay to lets then 0.01 per cent, end only the 54 win. beta 
activity remained. 

The majority of the data in this set of experiments wee 
obtained using a 100 millicurie radium-beryllium externel 
neutron source. This was the lergest source availeble at 


he time the experiments were begun. Five, 1 curie 
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plutonium-veryllium sources beoeme avalileble leter, and two 
sets of cata were obtained using this lerger neutron seurce. 
The first experimental errangenent consisted of the 6 in. 
lattice spseing with the 100 millicurie neutron source. A 
complete gurvey of the lower horizontal row of foils wlus the 
entire vertical row ef folls wam texen in this configuretion. 
Folis were svaced at 6 in. intervais from the east face to 
the centerline of the assembly. Next 8s verticel survey vas 
made using cadmium covered center line foils te determine 
tee cadmium ratio. The third set of data wes obtained with 
the seme Gin. lattice, but with only the £2 center line 
follse loaded. this wes to investigete the extent of flux 
Gepression ceused by the 152 folle of the orlginsl loading. 
The lattice sgpscing was then changed to 8.5 in. and 
with all foils loeded only the 18 in. and 350 in. (center Line} 
folls were counted. A set of folla wae also rleced commletely 
through the aseembly in the y direction. Dets obtsined from 
these folla were used in determining any flux variation due 
to the unsymsetrical foil loeding employed. 
Fuel loading wae changed again, and the 18 in. end 50 
in. flux measurements were repeated with a 12 in. lettice 
Spacing. 
Tre 1OU milliicurie radlum-beryllium neutron source was 
revleced by e Cluster of five, 1 curie plutonium-beryil ium 
neutron sources. lwo further eyvperimental runs were con- 


aueted. First a survey wes taken wit 911 fuel elements 
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removed. These Gata were used in Getermining the effeet of 
Subcritical muitipliestion. Finaliy, @ set of acta wes eol- 
lecteo usice the 6.5 in. lattice with the Larger neutron 
Source. These date provided s correiation tor the &.c& in. 
Lattice between charseteristics obtsined «ith the smell 
external neutron source end eherreteristics obtained with 

tee large eaternnl neutron souree. Again one row of Fells 
wes extended entirely throug® the sseembly to investigrte tre 
@yasetry in the y Alrection. 

i*e indium folis were positioned in e holder sei proxi- 
m@tely C.00 in. from the end window Gelger-j"ueller tube in a 
cunetant pesition reletive to the tube. Im area cese the 
araa Of the foll wae grester than t>e eren subtended by the 
enc winteow tube, and so correction wee mece for weleht vaeri- 
etione due to the alight area veriations among the foils. 

A atrontiugeyttrius esamle was uged a8 @ atamverd gouree, 
ena this semple wes counted st leest onee during exch ron. 
Thies yrocecure wes followed to insure thet no change in the 
experimental arrahgement or devietion in the Gelper Melier 
tube charseteriatics would creep in uncetected end invelidete 


the correlation between dsta cktreined on different ocecersions. 


















Oe ee ee! ot eee cee ee Cy : 
ee ie, ciel? —miodee ein 
ee ee i Vee ie bd me — 

‘et Re OA el my CETTE ies ate eee) 
OS 6 ID WD LO eee eed” 
ORD Gielige mobs ree ees Gr, es ee 
CAP i ee eee ie Pet | dans 
al Le ee et a ve) 4 
‘ eiawie © ee ap 
Ce Gales ©) Gail Wy cries <nome! 
8 et er lll le ee ee ce tee 
OF ee ie 6 8 Fy br leet Ls em 
8 Ce ce ee et ee, ee a ol 
AE ARETE ee et Le Om ke esl 
ee ee = ew od 
Ces a ed eb. @ 
— 2 abe cof a, a + ' 
ee es OT) Oe OD j 
Sitesorge om ce HY) * Hts — 
UMAR Ed Cee deededbame or por as Gn Ueeld ee 
ET RANG Oe Oe cee heey rere 














i 


= 


19 


DISCUSSION OF RESULTS 


fhe results of this series of flux measurements fre pre-~ 
sented graphically ae Figures 5 through 14. The deta on 
which these figures sre besed ere tebulated in the Apnendix. 
In a number of the figures s rorizentr] flux dietribution 
46 plotted. Since the theoreticel flux distribution hes the 
shepe of e cosine curve, cofxine curves were faired through 
the experimentel points rather than attempting to fit some 
other smooth curve te the pointe. 

& Getalled survey wes mace with the G in. lsettice and 
the gmell externel neutron source. Due to the emell activity 
of this source, counting rates were quite low end stenderd 
Geviations were necessarily large. Extended counting periods 
would have lowered deviations, cut this wes not voreetical 
Gue to the excessive ceeay of the 54 min. beta setivity while 
counting over s longer ceriod. Figure 3 disclays the flux 
Gistricution in the x direction as comprred to the theoretical 
ecvsine distribution. Vertical heig*t of the symbols indicetee 
g@tandsrd deviation. It is noted thet the seetter of experi- 
mental pointe decreases anc the flux dletribution conforms 
gore closely to the plotted cosine curves as y incresses 
toward the center line of the assembly. 

It was desired to determine whether a large number of 
indium folle present in the assembly produced a sersurable 


effect on the overall flux level. Figure 4 shows that any 
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verintion due to foil lording ie within the lisite of normel 
scattering, and only one cosine curve is plotted throug’ both 
pets of experimentei points. This cencinglion 16 slac sup 
ported by the deta of Figure &- Figure 6 illustrates the 
aymmetry in the y direction even thoug> seny core follsa yore 
pleeed in the exst half of the sssembly then in the weet half. 

In order to determine toe degree to which the faet source 
neutrons had been thermalilrec, a vertical flux survey was 
mace using cadmium covered Indiuw folls. Figure oO diselaye 
these dats contrested with the compereblie data obteined with 
bare indius follea. A constant reatto of these two fliuxes 
would indicate thet all of the reat souree neutrons >eé been 
thermaliced (4, r- 124). The date indieste thet feet neutrons 
6ti11 reaain fer some distenee uc from tee bottom of the 
@seembly. ib8 fpperently ehanging ratio neer tre top pey be 
ve to atatistical sesttering of the nolnta becruse of the 
very iow counting retes involved. It le believed trat further 
investigetion using the isarge external neutron source woulk 
enow this trend to be in error. 

Figure & eospares tre fliux obtained with each of three 


dattiace spacings. The &.o in. lattice gives the highest flux 


bevels of tre three lattices investigated. fome of the deta 


Obtelinea sith the ag in. inttice is believed to be in error, 
as these were the la'tt counts made before an equipment failure. 
Sluee the nature of the fellurea wes sueh thet the counting 


rate becawe extremely bigh just before comslete failure, 
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these experimental points ere believed to indicerte en un- 
peasonebly high flux. Therefore, the experimental points at 
xs £7 in. end greater were disregarded. Instead, the cosine 
curve wes fitted to points for x less then 27 in-, and the 
curve drawn symietricslly. The dashed line indicestes the 
extrepoleted portion of the curve. 

Frior to the completion of all experimentel ceterm' ne- 
tions a lerger neutron source became svelleble. Some ista 
were obteined using this larger source with tre 5.5 in. 
lattice. These results are presented both te confirm the 
findings obtained with the small neutron source end to demon- 
strate the great incresse in flux levels obtained by use of 
& larger neutron source. 

Figures 7 and & display the flux distribtution in the x 
ani y clrections, respectively, with eaeh of the neutron 
Sources in pleee. In both esses tre higner flux levels, 
corres: onding to the lerger neutron source, show little 
scatter of experisgental points, excellent symmetry, end a 
g00d conformance to the theoretical cosine distribution. As 
hes been noted previously, Figure 8 indicates no mersursble 
deviation from the cosine flux distribution due to unsynre- 
metrical indium foil leading in the y direetion. 

Qne set of cesta wes obteinesd with all fuel elements 
removed from the assembly. The flux levels present unter 


these conditions, when comprred with the flux levels present 
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in the 8.5 in. lattice, indicete the effect of subcritical 
multiplicetion. Figure 9 shows that within 35O in. vertically 
of the base of the essembly the empty grephite lettice hes a4 
higher flux level than Goes the one with fuel elements in 
place. his ig an indication thst more neutrons are being 
absorbed by the fuel than ere being produced by fission within 
the fuel. Above 30 in., however, the number of neutronea 
produced by fliesion in the urenium exceeds the number sbsorbed 
in the fuel, end the flux level with fuel elements in place is 
higher than tre flux level of the empty sasembly. Figure 10 
eomparee the flux level at the 50 in. pnlene for verlous 

velues of x horizontally. As steted stove, the overril 

effect of the fuel elements is essentielly zero et this dis- 
tance from the base. 

Une of the most importent chrracteristics which esn be 
determined experimentaily by meene of e subcritical assembly 
4a critical buckling. As developed previously, when flux (or 
&@ velue proportional to flux) 18 plotted on a logerithmic 
scale versus distance from the external neutroa source on a 
linear acele the Glepe of the line thus octained is ganma. 
if the physical dimensions of tre subcritical essembly ere 
Known, critical buckling cen then te determined by using 
Hyustien cl. 

In theory, eritical buekling cen be determined by first 
computing the infinite multiplicetion factor, kag, end then 


solving the transcendentsl equation 
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Figure 3. Effeet of suborition: multiplicrtion, 8.5 ineh 
Lattice, x = 27 inches, y = 50 inches 
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Figure 10. Effect of suberiticnal multiplicetion, ©.5 ineh 
lattice, »y «= 30 inches, 2 « 30 inches 


1 a L“B 
The infinite multiplicstion factor is computed from the 


four feetor formule 


Computing > from the relation 
om) = (24) 
> fuel 


gives 5 value fer natural urenium fuel of = 1.565. Cross 
sections used ere those given by Keplen (5, ». 4980). From 
&@ figure given by Classtone (6, ». 195), the fast fission 
factor, € , for the aize fuel elements used in the suberiticesl 
assembly is 1-050. ‘These two factors sre independent of 
latiice specing end remain the sme for e11 three latticec. 
Toe resonance eseare probability, 5, veries with the 
reletive amounts of fuel end moderator prevent, and it is 
therefore different for esch lattice specing. the method of 
Glagstone (6, pp- 191-192) was used to determine p. Velights, 
Volumes, and densities of fuel, cladding end graphite were 
Geterstined from representetive samples. The density of the 
graphite was computed to tre 1.56 g/en®. £8 overeil dimen- 
Sions of # fuel element were senagured to be 5.40 in. long by 
1-050 in. in diameter. Eecr fuel element weighed 2.018 kg, 
giving It en average Censity of 16.04 g/em”. Bach fuel 
element coiuteined a rod of naturel ursnium 8.00 in. long end 


1-00 in. in Giameter. This left eladéing of 0.040 in. thiek 


al 





Fs 


with G.200 in. end clates. Yf-e density of the cluminum 
eisdding end binder (eseumed the same) wes teken to be that 
ef sluminua, 2.71 fen. Thee density of the urenium was then 
Computed te be 19.0 g/cm. 

Using these dain, the squere Feterogeneour lattice wes 
treated of a aystem of equivelent evlindrice] unit cells (4, 
p- 265), and p wes detersined for each case. 

the Tourth feector, thermal utilizetion, f, else veries 


with Latuice specing- Thersial utilizstion wes esmputed os 


3 


shown by Murray (7, p- 37), ond apelin the equsre lettice was 
treated 18 a@ syste of equivelent cylindriecel unit cells. 
Allexance was wade for the presence of fuel, eledding, vir 
speee, and ypraphite moderatcr. for toe eslze fuel element 
used, the disadventsge factor, F, waa co@uted to be 1.089) . 
Computed values of m , €, b, FT, OHA Be, are tebnlated 
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teble 1. thtoretical calculstion of eritical buchiing 


























Value GS in- lettice 6.6 in. lattice lc in. inttice 
y 1-385 1.353 1-386 
é 1-050 1.030 1.056 
: 0.808 0.806 1.947 
f G.928 0.873 Q.9787 
ms 1.047 1 .0¢g 1.015 


B=(on~2} és x 1078 137 x 1078 8.5 x 10°" 
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daving Getermined k.., a first spproximrtion for eriticrl 
buckling can be obtained from Formula £5 which is v-1id for 


@ large reactor 


pB~ «kode k= 1 (25)(4, p- 216) 
we L*= ao 


Ferai ege, “(, for graphite is 360 om’. For a heterogeneous 


lattice, sage mey be taken a6 10 per cent grerter then that 
for the pure moderator (6, p. 201), giving as valwe of 585 cm’. 
Using velues of f just computed, a” wag obteinec from the 
formula 

L® 2 LE(a - r) (26)(4, mr. 260) 
Yalues of buckling, obtained by Equation 26 were then sub- 
etituted in the transcendental equetion, Equation 22. The 
meximum error wes 0.5 per cent. Slight adjustments geve 
Values of buckling correct within ©.1 per cent. 

The results of the above csleuletions for sill three 
lattices are tebuleted as Table 1. 

The regulte of experimental determination sre presented 
es Figures 11, 12, end 15. The bucklings computed from these 
figures sre listed in Teble 2. 

In the interest of clerity, #11 buckling dete hove been 
displayed on a common figure. Figure 14 ecomprres theoreticsl 
Duckling from Teble 1 with experimental values determined 
from Figures 11 through 15, as tebulsted in Teble 2. 

Figure ll compares vertical flux distribution for the 


three lsttices st x @# 27 ln. and y = 16 in. Ko atternt was 
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Figure ll. Lattice eotparison - vertienl, x « £7 inches, 
y = 16 inches 
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Plgure 12. Lattice comiprrison « verticel, x = “7 inches, 
y = 30 inehee 
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Figure 15. Vertical survey, &.5 ineh lattice, different 
neutron sources 
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Figure 14. Buckling versua insttice specing 
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Table <. Experimental values of ecriticel buckling 





Hemre a 
6 in. 8.5 in. 12 in. 
Conriguration lattice lattice lattice 





Small neutron source, ~ 6 Mot 
Xs 27 in-, y= 18 in. 141 x 107° 181 x 107” determined 


Smell neutron source, Ps o 
x #27 in., y= 30 in. 161 x 107° 65 x 10° is x 1076 


Large neutron source, Not Not 
x= 27 in., y = 50 in. determined 85 x 1076 determined 





made to estimate Buckling for the 12 in. lattice es the 
experimental pointe are too bedly eesttered. This is be- 
lieved due to the equipment failure nen tioned previously. 

Figure 12 preeente the same compsrisons et y = 30 in. 
This time the data for the 12 in. lattice apperred to be more 
reliable, and buckling was computed for #11 three lettices. 
Figure 15 presenta the data for the &.5 in. lettice from 
Pigure 12, and coupsres it with dste obtained for the came 
lattice with the large externel neutron source. Prrallel 
Streig¢ght linee have been ple tted for the two sets of data, 
end the experimental points fit theae lines very well. It is 
noted that the slope of theee lines gives s lower velue of 
buckling for the 8-5 in. lattice at y = 30 in. then did the 
Gete of Figure 11 taken at y = 18 in. 

It should be noted thet a very small chenge in the slope 


of a line makes s meorked difference in the velue of ortticensl 
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buexling computed from that slope. For this reason, reperted 
experiments at a high flux level would be necessary before 
completely reliable values of critical cuckling cculd be 
determined for a given lattice errangement. The experimental 
critical bucklings shown in Figure 14 should be considered to 
be indicative of general magnitude rather than highly securete 


velues. 
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CONCLUSIONS 


The opereting cheracteristicn of the euberiticel sseembly 
under investigation agree quite well with the secepted theory. 
Tre &-5 in. lattice is the optimam lsttice, resulting in 
higher flux levels end the largest velue of critics] buckling. 
All experimentel buckling determinstionsa proved to be posi-~ 
tive, which indlestes thet a full scsle resctor could be 
constructed of the same configureticn es either the 6 in., 

4.5 in., or 12 in. lattice. There 18 no espperent measurzble 
flux depression due to indium foll loading, even when foila 
are losded unsymmetricelly. Horizgontel flux €istributions 

ere symmetrical in totn the x end tree y directions, end within 
the limitatione imposed by staticties they sgree vell wit 

the cosine distributions predicted by theory. 

the euberiticel sssembly could be opersted setisfectorily 
using the smaller (100 millicurle) external neutron e#oures. 
However, the scattering of the experimentel points cbteined 
with this source 1s quite severe, end stetistical devintions 
ere very lerge. With the larger neutron source in plece, flux 
ijeveis cre higher by ao freetor of five or more, with corre- 
S@poncingly sieller devistions and less seettering of experi- 
Benteal points. 

The fast neutrons from the external source sare not 
entirely thermalized until they heave penetreted arproximately 


OO in. into the assembly. 
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This suberities] essembly, wit’ the lerge external neutron 
source in place, can be used to great advantnye as an iastruc~ 
tional tool. Although the experimentel results obteined were 
not extremely aceurate, it is believed thet the deta can be 
reprocuced with sufficient consistency snd sceurecy to cor 


relate eXperisenisal results with classroom caleulstions. 
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SUGGESTIONS FOR FURTHER STUDY 


Further study utilizing the seme suberiticel argembly 
ueead in this work might well be cerried out in one of three 
geherel directions. 

In the first of these categories lies the determination 
of various other properties of the asrembiy, such as diffu- 
Sion length, Fermi ege, end effective cross section. Methods 
ef cetersining these quantities experimentally sre given by 
Rurrey (7, p- 1G5). None of these properties has been deter- 
mined experimentally for thia perticuler suberitienl essembly. 

The second group of possible experiments incluées 411) 
those in which some besic physicel prremeter is veried end 
the effeet on nuclesr properties is observed. Exnoerimentsa 
euch a6 the study of local flux depressions under veriocus 
pertial fuel losds, simuletion of coolsnts, terpernture 
effeets on nuclesr prupertice, vaerisetion of neutron source 
geometry, and simuiction of econtrol rods all come under this 
Classification. Por sny lattice snaeing cther then 6 in. 
there is on appreciable volume of sir conteined in the emp ty 
fuel channels. Flux meseurements should tbe mede with these 
fuel channels fliled with grephite to determine the effect, 
if eny, of the empty channels. 

the leet type of experiment for which the suberiticel 
essembly might be used is one in which the eeserbly eets pas 


a testing ‘ecllity, rather then being investigeted itself. 
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This might inelude feaeibility studies of # new type of 
neutron counter, or the development of s moveble, continuously 
monitoring end recording survey instrument. This elso could 
include the experimental justificetion of newly ceveloned 
theory. An exemple of this type of applicrtion is the devel- 
opment of some new and leea critical method of experimentrlly 


detersining critical buckling. 
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Table 5. Key to tebles in Anrendix 








Neutron Lattice 
Table Run source spacing Comments 

4 au anell 6 in. 

S Ya enell 6G in. Cadmium covered 
6 ») emnil G in. Light foil load 
7 4 Small 8.8 in. 

3) s) emnell 12 in. 

B 6 large none All fuel resoved 
LO ? lerge &.5 in. 
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Tactle 4. (Continued) 





Saturated Sature ted 

Position activity(ce/im) Position activity(e/7) 
§-12-724 6 F=-6-18 50 
P-i2-30 6 Vm Bee 24 Fa 
Eel1-C CG B= 6=30 RS 
hL-L1l=-6 3 aha 14 
Belles KS Vi ww Bm & SP 
f-i1-18 10 We4-12 RQ? 
Rell-.4 & Fhe Ee 106 
£-11=50 8 bh ae 24 170 
Y-10-0 lL b whe 30 186 
aL Ge 4 P= 3-0 Li 
he-lvue-le 16 —: -§ 70 
&-10-18 thet” oe 117 
Selve-24 19 je 1 8 185 
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Saturated Jature ted 

fesition activity(c/m) Position ectivity(c/m) 
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Tevle 7. fiun 4 
Satureted Saturr ted 

Position activity(o/n) Position activity(c/m) 
B-Loelts 1 C-b—) 93 
B+ 13-50 3 G-5=-30 103 
b-12-15 y F-5=-18 105 
t-16.-30 Fd Fae b=—30 132 
h-11-185 6 Vw Se 1G 160 
h-lLie du 15 B= G~350 139 
b=-10-1% x E-4-18 of 
E~10—-36 15 Pe4—30 204 
h-9-12 ie C-3—18 ~ 18 
L=-9-30 4 Ew FeO 562 
B-G-16 39 TH-2=-18 210 
£-6—30 mH) T.-L =30 475 
b= 7-915 46 E-i-18 538 
pn 7 OD 65 =-1=-50 610 
1—-6-185 67 h.-8-0 12 
E-6=30 §3 Te G5 4 
fh-S-18 54 ie ey 74, 
fa~ 0-350 19 = Ga 4 1246 
b=6-16 46 Pe Be 36 129 
B-5-30 69 Fm Gah 2 120 
C- 5-165 7h Fe Hh -4 BO 
e~-§~30 115 Wn Be 54 48 
he Se 1S Os r~ 5-66 ye) 
5+ 8-30 1z9 

J-S-18 2s 

J-5=-50 18 

J-~5-15 56 

t=5-350 58 

w= 5-165 Via 

he b= 30 oF 








Yavle &. Hun 5 








cetureted 

Poai tion ac tivity(e/m) position 
Fe-2 5-18 3 JelJ5O 
E-jo-50 & T=5-18 
B-12-18 6 J~-5-5C 
oe wares) Es iw Se 1B 
LeL1-18 14 i= Se 36 
be-L41=-350 7 Go b=} 8 
h—-16-18 Li G-5=30 
N-i-50 17 FeS-1 8 
Vevei& 1& Pe $36 
Deve KL ae Fw Gm 1 A 
£-8-15 28 6-5-3090 
h-G-35. 46 hm S=-15 
ku 7-16 53 ra 4~ 50 
he 7= 50 49 Ye J- 1B 
EmO-1 8 6z Fe 5-30 
be Ge 3u 8 he l=15 
AnG-1 8 18 bei 30 
hwt— 30 LS E-1-18 
be SeL 49 Ve] 39 
Bev 350 5d 

CH Se-Lb 70 

C= 030 9O 

Lebel Ge 

tobe 50 1135 

J-l-18 14 








Setureted 
aetivityle/n) 


18 
Di. 
Sz 
9a 
102 


103 
1£% 
164 
145 
LLCO 


L3E 














Table G9. Run 6 
Saturated Setursted 

Position activity(c/nm) position eetivity(c/m) 
R-1L5-5u 9 Jo He 30 135 
Beli=-30 2 re G20) 70 
H-11-50 33 Fe5-8 235 
E=-10-50 O35 Fa fel 339 
f-9-50 73 Pe S18 555 
h=-5=50 149 Pa fe 4 650 
Eine Fm QO Lad Fe Be 50 AB4 
E=G=350 399 Te Ba 36 EG 
f-4-50 L107 Pe -42 587 
Re 5-50 1842 F- 5-49, 3E2 
be2=-30 2645 Kae Sie G4. LAO 
f=-]=350 5538 F-5-S0 76.5 
An 0-50 LZz9 

B= b= 50 316 

ee 469 

D~- v= 30 618 

Feu. 5 699 

G= 0-50 G17 

h=-GS= 30 479 

1-02-50 30 





Table 1. 


Run 7 








Seturated Satureted 

Position activity(c/m) fesition ectivity(e/m) 
f-1 3-50 £2 J-5=-350 141 
PelLi-3O 50 B- 5-6 85 
t-3)-350 61 te se) 244 
£~10=30 7B 4 =~ Ge} o S95 
B-9-350 110 P-5-1 8 53] 
F-BH.59 135 fim 5-4 BHO 
B- 7-30 ests) Fa 5-50 70 
F+6=39 432 Re Gan 5S S89 
£E-4-30 1094 Pa 5a 4 2 541 
Ee 5 30 16006 f— 5-48 394 
Ref =50 £330 Wem Fe H4 4 
Pm} = 50 ILLS - S= SO ES 
A=-6-30 136 

RW #~39 315 

C-5-30 492 

Aetert16) 5EE 

PuH.59 7% 

Ge b= 32 544 

fie be 30 406 

J -5-30 319 





Ce a I RD lp ci 
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